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I .  INTRODUCTION 

This  final  technical  summary  report  covers  the  period  of 
16  December  1965  to  1*J  December  1967.  The  success  which  we  have 
attained  in  the  realization  of  our  research  objectives  is  demon¬ 
strated  by  the  summary  of  aooompl ishments  under  this  grant  in 
Section  II.  Part  A  of  that  section  decorlhos  th®  experimental 
results  by  topic.  Part  B  is  a  list  of  publications  sponsored 
by  this  grant  (excluding  those  listed  in  the  1965  final  report). 
Part  C  is  a  list  of  talks  given  and  national  and  international 
meetings  attended  under  the  sponsorship  of  this  grant .  In  Sec¬ 
tion  III,  we  make  recommendations  for  future  work  in  this  field. 

II.  SUMMARY  OF  ACCOMPLISHMENTS  UNDER  THIS  GRANT  TO  DATE 

11965-1967) . 

Part  A.  Description  of  results  by  topic 

1.  The  main  accomplishment  of  this  program  to  date  is 
that  the  elastic  constants  and  their  pressure  and 
temperature  derivatives  have  been  measured  on  many 
specimens.  The  actual  measured  data  are  completed 
for  polycrystalline  AlgO^,  a  single  crystal  spinel 
(MgO*2.6  AlgOg),  single  crystal  garnet,  polycrystal¬ 
line  Mg2SiO^,  and  polycrystalline  ZnO.  This  is 
shown  in  Table  1.  In  addition,  the  data  on  poly¬ 
crystalline  CaO  and  polycrystalline  BeO  are  partial¬ 
ly  complete.  These  data  will  provide  future  refer¬ 
ences  for  many  thermodynamic  considerations  of 
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properties  in  the  earth's  interior. 

2.  The  experimental  approach  we  use  allows  us  to  mea¬ 
sure  the  value  of  the  bulk  modulus  K0,  and  the 
derivative  of  the  bulk  modulus  with  pressure,  K^. 
(Here  the  subscript  zero  means  evaluated  at  zero 
pressure.)  We  have  found  that  to  a  high  degree  of 
precision,  K'  is  independent  of  pressure  up  to  the 
limits  of  our  experiments,  3  kb.  We  also  find  that 
if  we  assume  that  K'  is  independent  of  pressure  at 
much  higher  pressures,  we  can  predict  the  results 
of  density  pressure  measurements  obtained  by  other 
kinds  of  experiments.  The  basic  for  predicting  ex¬ 
perimental  values  Js  that  since  K  is  a  differential 
of  pressure  with  density,  the  assumption  is  a 
constant  yields  a  differential  equation  which  in¬ 
tegrates  to  a  particular  p-P  relationship.  This  is 
the  Murnaghan  equation  of  state,  which  contains  only 
the  two  measured  quantities  K0  and  K^,  ard  no  arbi¬ 
trary  constants.  As  a  consequence,  our  acoustic 
measurements  have  provided  a  method  of  tying  into 
and  calibrating  shock-wave  data  measurements  at 
very  high  pressures.  This  method  of  correlating 
shock  data  and  acoustic  data  was  not  used  before 
our  work  began.  We  therefore  believe  our  work  has 
had  an  impact  on  shock-wave  analysis. 
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3.  By  using  the  equation  of  state  resulting  from  the 
condition  «  constant,  we  have  been  able  to  de¬ 
rive  equations  for  the  seismic  velocity  as  a  func¬ 
tion  of  pressure,  and  the  coefficient  of  thex'mal 
expansion  as  a  function  of  pressure.  Because  of 
arguments  discussed  in  Part  2  abovo ,  we  to]1eve 
that  these  equations  are  valid  over  px-essm-c  rnn^os 
of  several  hundred  kilobars  (in  the  absence  of  phase 
changes).  These  equations  are 


density, 


seismic  parameter. 


(Ref.  2,  p.  7) 


(Ref.  2,  p.  7) 


thermal  expansion,  0  ^  a0  ^1  +  Kq  (Ref.  10,  p.  8) 

The  important  point  about  these  equations  is  that 
the  parameters  pQ,  aQ,  KQ ,  and  are  all  determined 
at  very  low  pressures,  so  that  values  of  constants, 
such  as  reported  in  Table  1,  are  sufficient  to  esti¬ 
mate  the  quantities  p,  $,  and  a  at  very  high  pressures. 
These  equations  are  useful  to  theories  concerning  the 
history  and  structure  of  planets.  No  claim  is  made 
for  their  being  formally  rigorous,  but  the  equation 
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for  p/  seems  to  be  as  accurate  as  the  best  avail- 
Po 

able  measured  data;  and  if  that  equation  is  valid, 
the  two  equations  which  follow  are  also  valid. 

4.  We  have  found  extrapolation  formulas  for  the  elastic 
constants  at  very  high  temperatures.  (This  original 
work  was  not  credited  to  this  grant  because  it  was 
done  with  equipment  at  Bell  Telephone  laboratories 
while  Dr.  Soga  was  on  leave  from  Lamont  to  Bell  Labs. 
It  has  become  very  important  to  this  grant.)  Ac¬ 
cording  to  our  theory,  as  confirmed  by  experiments 

up  to  1300°C,  the  variation  of  bulk  modulus  with  tem¬ 
perature  at  high  temperatures  is  given  by 

dK  «yCp  (Ref.  5,  p.  7) 

If  T —  (Ref.  7,  P.  8) 

where  6  and  y  are  measurable  constants,  V  is  the 
specific  volume,  and  Cp  is  specific  heat.  This  ex¬ 
pression  predicts  that  dK/dT  does  not  change  very 
much  with  composition.  This  is  demonstrated  by  mea¬ 
sured  values  in  Table  1.  An  important  result  of  this 
work  is  that  it  allows  one  to  predict  the  behavior 
of  velocity  with  temperature  in  the  mantle. 

5.  We  can  predict  the  variation  of  (3T/3P)  at  constant 
♦  -  /K7o  as  a  function  of  pressure  and  temperature 
at  high  pressures  and  temperatures.  This  results 
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from  the  extrapolation  procedures  outlined  in  (4)  and 
(5)  above.  This  parameter  is  very  important  in  the 
estimates  of  thermal  gradients  in  the  mantle  mater¬ 
ials.  Ac  an  example,  the  extrapolated  values  of 
( 3T/3P)  at  constant  are  given  below  in  degrees/ 
kilobar  for  the  ease  of  a  garnet. 


p 

T 

25°C 

500°c 

1000°C 

1500°C 

0  kb 

27 

21 

19 

18 

30  kb 

26 

20 

18 

17 

50  kb 

19 

17 

16 

100  kb 

23 

18 

16 

15 

We  are 

able 

to  predict 

the  variation  of 

the  critical 

geothermal  gradient  (the  gradient  at  which  <t>  is  con¬ 
stant  with  depth)  on  a  given  P-T  profile.  The  P-T 
range  available  for  this  extrapolation  extends  to 
the  limits  of  the  upper  mantle  of  the  earth. 

Part  B.  List  of  publications  (excluding  those  listed  in 
the  1965  final  report)  sponsored  by  this  grant. 

PUBLISHED; 

1.  "Temperature  dependence  of  the  velocity  derivatives 

of  periclase,"  J.  Geophys.  Res.,  71 ,  3007-3012, 
1966,  by  E.  Schreiber  and  0.  L.  Anderson. 

2.  "Seismic  parameter  :  Computation  at  very  high 

pressure  from  laboratory  data,"  Bull.  Selsmol. 

Soc .  Am. ,  £6/3)  >  725-731*  1966,  by  0.  L.  Anderson. 
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3.  "A  proposed  law  of  corresponding  states  for  Oxide 

compounds, "  j.  jQeophys.  Res.,  71,  4963-4971,  1966, 
by  0.  L.  Anderson. 

4.  "Pressure  derivatives  of  the  sound  velocity  of  poly¬ 

crystalline  forsterite  with  6%  porosity,"  J .  Geo- 
phys.  Res..  72,  760-764,  1967,  by  E.  Schreiber  and 
0.  L.  Anderson. 

5.  "Derivation  of  Vachtman's  equation  for  the  temper¬ 

ature  dependence  of  elastic  moduli  of  oxide  com¬ 
pounds,"  Phys.  Rev.,  144(2),  553-557,  1966,  by 
0.  L.  Anderson  (submitted  under  Bell  Labs  byline). 

6.  "Pressure  derivatives  of  elastic  constants  of  single 

crystal  MgO  at  23°C  and  -195- 8°C,"  J.  Am.  Ceram. 

Soc . ,  49,  404-409,  1966,  by  Orson  L.  Anderson 
(submitted  under  Bell  Labs  byline). 

7.  "High  temperature  elastic  properties  of  polycrystal¬ 

line  MgO  and  AlgO^,"  J.  Am.  Ceram.  Soc.,  49 .  355- 
359,  1966,  by  N.  Soga  and  0.  L.  Anderson  (submitted 
under  Bell  Labs  byline). 

8.  "Variable  air  transformer  for  impedance  matching," 

Rev.  Scl.  Instr.,  37 .  1625-1626,  by  P.  Mattaboni 


and  E.  Schreiber. 
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IN  PRESS; 

9.  "Elastic  moduli  of  single  crystal  spinel  at  25°C 

and  to  2  kbar,"  Am.  Ceram.  Soc .  .  by  F.  CcLrdher. 

10.  "An  equation  for  thermal  expansivity  in  planetary 

interiors,"  J.  Geophys.  Res.,  by  0.  L.  Anderson. 

11.  "Sound  velocity  in  rocks  and  minerals,"  chapter  in 

a  book.  Physical  Acoustics,  IVB,  Academic  Trees, 

0.  L.  Anderson  and  R.  C.  Liebermann. 

12.  "Elastic  constants  of  oxides  used  to  estimate  the 

properties  of  the  Earth’s  interior,"  Proceedings — 
NATO  Advanced  Study  Institute:  on  the  Application 
of  Modern  Physics  to  the  Earth  and  Planetary  .Inter¬ 
iors^  29  March-4  April,  1967,  Newcastle  upon  Tyne, 
England,  by  0.  L.  Anderson  and  R.  C.  Liebermann. 

Part  C .  List  of  talks  sponsored  by  this  grant 

Presentation  of  contributed  papers  before  the  1966 
and  1967  April  meetings  of  the  American  Geophysical 
Union.  Invited  seminars  at  Yale  University,  M.I.T., 
U.C.L.A.,  University  of  Hawaii,  Rice  University, 
California  Institute  of  Technology,  and  Air  Force 
Institute  of  Technology,  Dayton,  Ohio,  and  AFOSR, 
Washington,  D.  C. 

Invited  paper  at  NATO  Conference  on  Application  of 
Modern  Physics  to  the  Earth  and  Planetary  Interiors, 
29  March  1967*  Newcastle. 


f 


8. 


Invited  paper  at  fall  meeting  of  the  Daolc  Science 
Division  of  the  American  Ceramic  Society,  Pitts¬ 
burgh,  Pa.,  October  1965. 

TIT.  Recommendations  for  future  work 

The  results  described  in  the  previous  section  indicate  the 
desirability  and  the  rewards  of  continuing  our  work  and  extending 
it  to  a  new  set  of  objectives. 

OBJECTIVES 

Part  A.  Continue  the  work  of  the  previous  grant  (see  ob¬ 
jectives  listed  in  proposals  of  1$  December  l^^ 
and  15  December  1966)'. 

1.  Prepare  and  measure  the  sound  velocities  as  a  function 
of  pressure  and  temperature  on  unmeasured  oxide  com¬ 
pounds  and  minerals:  in  particular,  those  oxides  con¬ 
taining  transition  elements  Pe,  Mn,  and  Ni. 

2.  Estimate  the  limits  on  geothermal  gradients  in  the 
earth's  mantle  from  new  and  old  data  obtained  from 
work  on  this  grant. 

3.  Choose  the  experiments  and  interpret  the  results  to 
bear  upon  current  geophysical  problems  such  as  earth¬ 
quake  mechanisms,  convection  in  the  mantle,  and  sea¬ 
floor  spreading. 

l*.  Find  the  thermodynamic  variables  in  the  neighborhood 
of  phase  changes. 

Part  B .  Begin  new  work  on  the  measurement  of  thermal  ex¬ 
pansivity  as  a  function  of  pressure  and  temperate 

1.  Measure  the  thermal  expansivity  in  oxide  compounds  a 
a  function  of  pressure  at  different  temperatures. 
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Part  C.  Find  practical  wavs  of  measuring Uie__elastlc_con- 
stants  and  their  pressure  de r 1 v a 1 1 ves  of  the 
monocllnlc  crystal  class 

Part  D.  Exploit  the  spherical  resonance  technique  for 
measuring  elastic  constants 

1.  Compute  the  expected  behavior  for  spherical  reson¬ 
ance  for  spheres  made  from  nonisotroplc  substance. 

2.  Measure  the  shear  velocity  of  small  spheres  as  a 
function  of  temperature  over  wide  temperatures. 

Part  E .  Find  suitable  rock  specimens  for  laboratory  mea¬ 
surements  which  have  come  to  the  surface  from 
deep  In  the  earth's  Interior 

Items  1,  2,  and  3  of  Part  A  above  are  extensions  of  our 
previous  work.  Part  A*!  has  not  been  indicated  before,  so  some 
brief  comments  will  now  be  made.  Suppose  there  is  a  phase  change 
at  pressure  P*.  Acoustic  measurements  will  predict  the  volume  at 
pressures  above  P#  by  the  extrapolation  techniques  described.  By 
using  shock-wave  techniques  the  volume  above  P*  can  be  measured, 
and  can  be  interpolated  below  P#  to  zero  pressure.  Thus  one  can 
estimate  the  volume  change  aV  between  the  two  phases  above  and  be 
low  the  transition  pressure.  The  product  of  the  pressure  and  the 
AV  is  the  change  in  enthalpy  between  the  two  phases.  By  using 
the  condition  that  the  free  energy  difference  vanishes  at  the 
transition,  and  by  exploring  new  (yet  undefined)  techniques  for 
measuring  the  entropy  difference  at  all  pressures,  one  may  hope 
to  arrive  at  the  basic  information  for  thermodynamic  constants 
relative  to  phase  changes.  This  point  should  be  explored  theo¬ 
retically  and  experimentally.  Such  information  may  possibly  lea 
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to  the  prediction  of  phase  changes. 

Part  B  of  OBJECTIVES  concerns  thermal  expansivity.  In  the 
course  of  our  work,  we  have  found  that  the  data  on  thermal  ex¬ 
pansivity  are  important  to  transform  and  interpret  our  acoustic 
data.  V’e  have  measured  our  own  thermal  expansivity  data  as  a 
function  of  temperature.  Very  little  work  has  been  done  since  the 
time  of  Bridgman  on  the  measurement  of  thermal  expansivity  as  a 
function  of  pressure.  This  work  should  begin.  The  relation  be¬ 
tween  acoustic  data  and  the  a-P  relationship  is  given  by  the  thirc 
equation  on  page  3. 

Part  C  of  OBJECTIVES  concerns  the  problem  of  making  acoustic 
measurements  on  samples  of  low  symmetry.  The  number  of  data  whic) 
one  must  measure,  on  the  basis  of  present  theory,  is  formidable. 
The  theoretical  aspects  of  this  problem  should  be  considered  in 
order  to  make  an  experiment  program  practical. 

Part  D  of  OBJECTIVES  concerns  the  spherical  resonance  tech¬ 
nique  which  has  been  exploited  by  Soga  and  Anderson  (J.  Oeophys. 
Res. ,  72 ,  1733-1739,  1967)  as  a  method  for  measuring  small  tek- 
tites  on  a  NASA  project.  In  the  future,  it  may  be  useful  for 
sound  velocity  measurements  on  small  but  rare  minerals.  It  thus 
may  be  pertinent  for  the  Apollo  lunar  material  program,  if  that 
program  ever  materializes.  However,  we  are  concerned  here  with 
the  use  of  spherical  resonance  to  general  problems  in  geophysics 
and  rock  mechanisms.  One  advantage  of  this  method  is  that  sphe 
can  be  prepared  easily j  while  the  main  disadvantage  is  that  the 
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theory  has  not  been  worked  out  for  spheres  of  arbitrary  crystal 
symmetry.  We  want  to  solve  this  important  problem,  so  that  we 
may  be  provided  with  a  very  powerful  technique  for  dealing  with 
small  specimens.  Lower  symmetry  causes  mode-splitting,  but  so 
does  aspherlcity  of  the  specimen.  These  lead  to  departures  from 
the  isotropic  case  and  must  be  distinguished  from  each  other. 

The  application  pertinent  to  our  work  is  the  adaptation  of  the 
spherical  resonance  technique  to  determining  the  effects  of  tem¬ 
perature  on  sound  velocity.  It  may  also  be  feasible  for  use  with 
pressure  experiments  as  torsional  mode  should  not  be  damped  by 
fluid  at  high  pressure.  Thus  we  would  expect  to  be  able  to  mea¬ 
sure  the  effect  of  pressure  on  the  shear  velocity  of  a  small 
specimen. 

Part  E  of  OBJECTIVES  concerns  the  acquiring  of  real  materia] 
from  the  deep  parts  of  the  earth.  Our  program  has,  from  the  firs 
been  designed  to  start  from  the  simplest  oxide  compounds  and  grac 
ually  work  to  more  realistic  components  of  the  earth's  interior. 
We  are  about  ready  to  begin  measurements  on  natural  rocks  such  a i 
eclogite,  dunite,  peridotite,  pyroxinite,  and  their  constituent 
minerals.  The  preliminary  survey  of  work  on  these  compounds  is 
tabulated  and  discussed  in  a  VESIAC  State-of-the-Art  Report, 
"Sound  Velocities  in  Rocks  and  Minerals,"  (by  0.  L.  Anderson  and 
R.  C.  Liebermann,  issued  November  1966  by  the  University  of  Mich 
igan  Geophysics  Laboratory  for  ARPA).  V/e  are  presently  looking 
for  suitable  specimens  for  our  existing  techniques. 
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There  are  two  field  problems.  First,  we  want  to  find  the 
available  rocks  from  the  deepest  regions  of  the  earth  -end  for 
this  we  need  to  apply  the  techniques  of  geology;  second,  we  want 
to  select  from  these  rocks  those  which  are  suitable  for  measure¬ 
ments  by  our  techniques.  This  involves  field  trips,  preferably 
with  the  help  of  a  trained  field  geologist. 

In  the  past  eighteen  months,  we  have  gone  on  two  field  trips 
to  find  specimens  using  funds  from  this  grant.  We  found  an 
aphanitic  nonvesicular  basalt  near  Lava  Falls,  Grand  Canyon.  We 
We  found  some  eclogites  at  the  Mule  Ear  kimberlite  pipe  in  south¬ 
ern  Utah,  which  apparently  were  derived  from  great  depth.  This 
encourages  us  to  seek  our  own  specimens.  We  propose  to  look  for 
suitable  samples  from  three  possible  kinds  of  locations.  First, 
xenoliths  from  volcanoes  (such  as  found  on  Mauna  Loa,  Hawaii); 
second,  ultrabasic  rocks  from  orogenic  zones  (such  as  those  occur 
ring  in  the  Franciscan  formation  of  California),  and  third,  ultra 
basic  rocks  from  cratonic  zones  (such  as  the  Mule  Ear  diatreme, 
Utah. 
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Orson  L  Anderson 
Principe  I  Investigator 
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